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Abstract 

O ' We propose a universal description of dark energy and modified gravity that includes 

all single-field models. By extending a formalism previously applied to inflation, we con- 
sider the metric universally coupled to matter fields and we write in terms of it the most 
general unitary gauge action consistent with the residual unbroken symmetries of spatial 
diffeomorphisms. Our action is particularly suited for cosmological perturbation theory: 
the background evolution depends on only three operators. All other operators start at 
least at quadratic order in the perturbations and their effects can be studied indepen- 
dently and systematically. In particular, we focus on the properties of a few operators 
which appear in non-minimally coupled scalar-tensor gravity and galileon theories. In 
this context, we study the mixing between gravity and the scalar degree of freedom. We 
assess the quantum and classical stability, derive the speed of sound of fluctuations and 
the renormalization of the Newton constant. The scalar can always be de- mixed from 
gravity at quadratic order in the perturbations, but not necessarily through a conformal 
rescaling of the metric. We show how to express covariant field-operators in our formal- 
ism and give several explicit examples of dark energy and modified gravity models in 
our language. Finally, we discuss the relation with the covariant EFT methods recently 
appeared in the literature. 
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1 Introduction 

The discovery of cosmic acceleration has motivated a tremendous amount of observational 
and theoretical activity. On the observational side, a wealth of precious cosmological data is 
awaiting us, thanks to the next generation of large scale structure surveys such as EUCLID [TJ2] 
and BigBOSS [3J. By accessing a very large number of modes, such experiments will be 
sensitive to the details of structure formation and therefore offer the unique possibility of 
discriminating between competing models of dark energy and modified gravity (hereafter, in 



2 



short, DE) on the basis of dynamical and clustering properties jl]. On the theoretical side, 
the number and variety of proposed models for the cosmic acceleration (see e.g. [SJE] for 
recent reviews) is by no means less impressive. Such an overwhelming production is certainly 
a resource for cosmology, but also a somewhat embarrassing load of material difficult to deal 
with. One reason for concern is the unnatural complexity of many DE models — especially 
when it is a naturalness problem that they are allegedly addressing. This is of course, to 
some extent, a subjective issue, difficult to quantify and find general agreement on. If we 
decide to maintain a democratic view on all different proposals and leave the final verdict to 
observations we face a more pragmatic problem: efficiently discriminating between different 
models. Each one has a certain number (and type) of parameters and has to be fitted against 
data independently from the others, by first solving the background and then the perturbation 
equations. As a matter of fact, in most cases, the performances of each DE model are compared 
"one against one" with those of the best contender so far, ACDM. 

This situation calls for a "unifying" and "effective" approach to DE modeling. "Unifying" 
in the sense that it should incorporate as many different models as possible as special cases; 
"effective" , here, in the non-technical sense of being readily testable by observations and not 
committed to specific DE models, nor to their original motivations. A substantial effort has 
been put in the recent years in this direction; see for instance [THIS]. 

Most models of dark energy and modified gravity can be described — in their relevant re- 
gimes — with the only addition of a scalar degree of freedom to the Einstein-Hilbert action. 
This is reminiscent of inflation, where a scalar field is needed to break de Sitter invariance. In 
the context of inflation, the goals mentioned above were efficiently achieved by the 'effective 
field theory (EFT) of inflation', initiated by Creminelli et al. in Ref. [T7] and then more 
systematically developed by Cheung et al. in Ref. |18^l . The idea is to apply EFT directly 
to cosmological perturbations, by treating them as the Goldstone bosons of spontaneously 
broken time-translations. The often invoked analogy is with the spontaneous breaking of 
the SU{2) x U{1) gauge symmetry in the Standard Model. In unitary gauge the would-be 
Goldstone bosons are "eaten" by the longitudinal degrees of freedom of the vector bosons W± 
and Z. At the price of loosing manifest gauge invariance, one can deal directly, at the EFT 
level, with the observable (and now finally all observed) low-energy degrees of freedom of the 
theory: three massive vector bosons and one Higgs particle. 

The use of the unitary gauge brings into cosmology similar advantages. The operators can 
be organized straightforwardly in powers of the number of perturbations; as a consequence, 
operators beyond the linear order do not affect the background evolution. Moreover, the terms 
in the expansion have direct observable implications. For instance, the cubic operators of the 
EFT of inflation formalism can be straightforwardly related to the observable three-point 
functions of the CMB HBJE2H23]. 

The application of this formalism to dark energy was initiated in [9], where generic minimally- 
coupled single-field models were considered. The main difference with inflation is that in the 
late time universe matter species are present. In this paper we extend the study of [H] to a 
non-minimally coupled scalar field in order to develop a unifying theoretical framework aimed 
to include all dark energy and modified gravity models. The logic of our approach can be 
summarized in the following two basic steps: 

1 This formalism made its first appearance in Ref. [T5], where it was used to study the coupling of the ghost 
condensate with gravity. See also for later generalizations to multi-field inflation and dissipative effects 

in inflation. 
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a) We assume the validity of the weak equivalence principle (WEPjl and therefore the exis- 
tence of a metric g^ u universally coupled to matter fields through an action S m \g ia ,, ip, 



mj i 



b) We write the unitary gauge action, i.e. the most general gravitational action for such a 
metric compatible with the residual symmetries of unbroken spatial diffeomorphisms. 



As emphasized in [171118] . the last point allows, although does not postulate from the beginning, 
the presence of a scalar field <j) in the DE sector. In unitary gauge such a scalar does not appear 
explicitly because it is "eaten" by the metric. More practically, the time coordinate is chosen 
such that it is a function of the scalar field, so that the fluctuations of <fi around the background 
vanish, S(f>(t,x) = (fi(t,x) — 4>o(t) = 0. The scalar degree of freedom can reappear explicitly 
in the action, together with full diffeomorphism invariance, by performing the "Stueckelberg 
trick", i.e. by performing an infinitesimal time diffeomorphism t — >■ t — ir(x) (see Sec. (12. ip ). 
where ir is now the field perturbation encoding the scalar dynamics of DE. 

In order to show the unifying power of this approach, let us write the unitary gauge DE 
action in terms of the universally coupled Jordan metric g^ u . This reads 



S = d A x 



M. 
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-f(t)R - A(t) - c(t)g 00 



+ S { $ E , (1) 



where R is the Ricci scalar, /, A and c are functions of the time coordinate t, g is the upper 

(2) 

time-time component of the Jordan frame metric and S DE indicates terms that start explicitly 
quadratic in the perturbations and therefore do not affect the background. More explicitly, 
the last three lines of eq. (j3J) below. In this paper, M* is the "bare" Planck mass. Its relation 
with the observed gravitational constant is given in some specific cases in eqs. (1541) . fl65l) and 

fl7DD. 

Let us first consider the square brackets. The time-dependent coefficient f(t) in front of 
the Ricci scalar is reminiscent of the non-minimal coupling function in scalar-tensor theories, 
except that here it is written in unitary gauge. The presence of this function is the main 
difference with the EFT of quintessence approach [9]. In the absence of matter fields, as 
discussed in [TS], f(t) can be re-absorbed by a conformal transformation of g^. Here, however, 
the matter sector uniquely singles out the Jordan metric g^, in which test particles follow 
geodesies. In this paper we mostly privilege the Jordan frame but we provide the reader with 
the equivalent, and often more handy, "Einstein" -frame version of our formalism in Sec. [5j 

The other two functions of time, A(t) and c(t), have also a very intuitive origin in the 
simple case of a canonical scalar field. For instance, in unitary gauge the kinetic scalar term 
{d<p) 2 becomes 

- l -{d<pf = —fTd^ -> -c(t)g 00 , (2) 

with c(t) = <Pq/2. However, the operator — c(t)g 00 can also receive contributions by higher- 
order operators such as (<90) 4 . In general, A(t) and c(t) can be expressed in terms of back- 
ground quantities such as the DE background energy density and pressure, respectively pr> 



2 With this choice, we aim to emphasize the "modified gravity" content of the theory and disentangle our 
analysis from the interesting but somewhat off-topic rich phenomenology of WEP violations. The choice of a 
universally coupled metric is stable under radiative corrections in the matter sector [25H27] , meaning that WEP 
violations are generally expected to be delivered by Planck suppressed operators. See Sec. [5] for a relaxation 
of this assumption. 
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and pd-i and the Hubble rate H. In Sec. I2.2[ by using the Einstein equations we find 



1 



2 



1 



A(t) = - PD-PD + M^(5Hf + f) , c(t) = - p D +p D + Mi{Hf-f) . (3) 
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The background dark energy density and pressure are defined by the modified Friedmann 
equations (l2Tj) and (122]) below and are thus derivable from the observed DE equation of state. 

Of course, structure formation will be also sensitive to the quadratic and higher-order terms 
contained in S EE and these will vary from one DE model to the other. In our EFT approach, 
the building blocks by which these operators are constructed are perturbations of quantities 
which are invariant under spatial diffeomorphisms, such as the upper time-time component of 
the metric g 00 , the extrinsic curvature of uniform-time hypersurfaces K^ v (esplicitely defined in 
eq. (jfJJ)) and the Riemann tensor R^ap [18]. We are also allowed to contract these quantities 
with the metric or among each other to form scalars and to take their derivatives. These 
perturbations appear in quadratic or higher-order combinations so that they do not affect 
the background and can be studied — and fitted against observations — independently and in a 
systematic way. Despite the richness of S DE , the higher-order operators are classified according 
to their perturbative effects so that at the linear order in the perturbations only a finite number 
of them have to be considered. We discuss all these higher-order operators in details in Sec. 12.31 
and below. 

In summary, including higher-order terms our Jordan frame action reads 



oo 



S = \ j d A x^g- [m 2 JR - Pd+Pd - M 2 (5Hf + /) - [p D + p D + M 2 (Hf - /)) g 

+ M*(5g 00 ) 2 - m\ 5g 00 6K - M 2 5K 2 - M 2 SK^SK^ + m 2 2 h^d^g 00 d u g 00 
+ X 1 5R 2 + \ 2 5R^5R^ + n\5g m 5R + nC^C^ + i^C^ kX C^ kX 

+ ^(Sg o y-mU6g o y6K + ... 

where 5g 00 = g 00 + 1, SK^ is the perturbation of the extrinsic curvature, SK its trace, 5R 
and 5R^ U the perturbations of the Ricci scalar and tensor, respectively, and C^ upa is the 
Weyl tensor. Mi, rrii, Mi, rhi and /ij are mass parameters while Aj and 7j are dimensionless 
parameters; they can all depend on time. 

Here follow some explanations, comments and a summary of results. 

• Set up. The action (j4j) is especially set up and arranged for cosmological perturbations 
on a FRW background. Every term can be expanded to arbitrary order. However, all DE 
operators contributing up to linear order in the perturbations are contained in the first line 
of eq. (jlj) and their coefficients are fixed by background quantities only, p D , p D , H and 
/, which can be fitted against observations (Sec. I2.2p . Quadratic terms (second and third 
lines) and higher (forth line, etc.) leave the background unchanged. Their time-dependent 
coefficients, Mf, fh\ etc., are sensitive to tests which depend on the clustering or screening 
properties of DE. 

• Jordan frame. The action (jl]) is written in the Jordan frame g^ v that universally couples 
to the matter fields. We privilege the Jordan frame for two reasons: it is more directly 
related to observations [27] and it is univocally defined by the coupling to matter, once we 
postulate the WEP. 
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• The scalar dynamics and mixing. The action (j3J) is written in unitary gauge and therefore 
only depends on the metric and its derivatives. Full diffeomorphism invariance can be 
restored by introducing a scalar degree of freedom ir via the Stueckelberg trick (see, however, 
Sec. I3.ip . after which any other gauge can be chosen. At sufficiently high energy, only 
the kinetic terms are important in the quadratic action. In Ref. [18] this is called the 
decoupling limit, because in the simplest CcLSGS clS well as in the notable one of Electroweak 
SU(2) x U(l) theory — the scalar decouples from the gauge fields. More generally, in Jordan 
frame the scalar can be kinetically mixed with gravity, in which case there is no decoupling 
limit. In Sec. [3] we study the mixing due to / and m\ in Newtonian gauge. We show that 
one can de-mix it and gravity by an appropriate local field redefinition, at the expense of 
coupling 7T to matter. 

• Comparison with EFT of inflation. We are considering few more operators than the ones 
included in [18] for the inflationary case. In the absence of matter fields, some of them (/, 
Ai, A2) can be reabsorbed by a field redefinition of the metric tensor. Others (like m\ and 
Hi) were not explicitly included and generally implied by ellipsis; we have included them 
here because we found them essential for some DE models. As mentioned, we also gave more 
relevance to the operator ml (there called Mf and also discussed in [9]), as it is responsible 
for kinetic mixing. 

• Universality. The action (jlj) describes virtually all DE models based on a single degree of 
freedom, as we detail in Sec. H] with several examples. What this action does not describe 
are genuinely higher-dimensional regimes of, say, brane scenarios, that cannot be encoded in 
a 4-d description, massive gravity away from the decoupling limit, theories including vector 
fields participating to the gravitational sector such as TeVeS (see e.g. [5] for references), as 
well as modifications of gravity of the geometrical type, such as those suggested in [291 [3D] . 
As discussed in Sec. 15. 1[ special and different couplings to the dark matter sector (like those 
needed for models of "coupled quintessence", e.g. [3TH33] ) can also be included with minor 
modifications. 

• Action principle vs. linear equations and regime of applicability. Despite the technical 
complications related to the gravitational sector (e.g. having to solve the constraints [28]) 
the formulation in terms of an action is completely general and therefore preferable to the 
linearized equations. The regimes of applicability of (jlj) span from the weakly coupled 
largest cosmological scales down to the UV cutoff of the EFT. This can be estimated by 
computing the energy at which unitarity is violated [TBIIM] . 

Other proposals employing an EFT approach to cosmic acceleration have been put forward 
in the literature, such as those of Refs. [TU ] [T2" ] |3"5] . They are based on a covariant description of 
a canonical scalar field together with higher-order terms in a covariant expansion up to four 
derivatives, in the spirit of the EFT of inflation a la Weinberg [36] . This approach is certainly 
a promising starting point in the direction of unifying models with an effective description. 
Moreover, the familiar language of (covariant) EFT, with the relative importances of mass 
scales and operators clearly displayed, can be used here in all its power to explore the space of 
DE models. As shown in this references, the model space can be tighten up by perturbative 
field redefinitions: models that look different from each other may in fact be equivalent EFTs 
at a closer scrutiny. 
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However, we would like to stress the advantages of our approach as compared to a the 
aforementioned one, or more generally as compared to the use of a covariant description using 
the language of 0. Most of these points are also well explained in [9"|rT5]l!2"3"]. 

• Covariance vs background/perturbations splitting. In the <fi language, any further term 
in the EFT expansion corrects the FRW background. This means that including a new 
perturbative correction boils down, in practice, to studying a new model altogether ab initio: 
solving new equations for the background, studying perturbations around it in order to seek 
for specific dynamical effects, etc. This is avoided in our approach, where perturbations are 
expanded about a fixed background once for all. Note however that this does not mean, 
as sometimes stated, that our approach cannot address the background dynamics: given 
action fl4]) with a finite number of operators, we can always construct a covariant action in 
the language of <fi which correctly reproduces the evolution of the background and of the 
perturbations, as shown in [5] and in Sec. 12.41 

• Stability. The stability and consistency of any model can only be evaluated once one turns to 
and g^y fluctuations expanded around a FRW background. This is also the case of theories 
that are known to have second-order field equations, such as the Horndeski theory [371 - 139] . 
but whose stability depends on specific signs of the quadratic kinetic Lagrangian of the 
fluctuations. 

• Relevance of operators. The relative importance of the various operators can be assessed 
only after expanding the action around a background. For instance, there are cases in 
which covariant operators that are naively higher order in perturbation theory become in 
fact relevant and should not be treated as perturbations in the standard EFT sense. For 
example, in the ghost condensate mechanism [19] a vanishing speed of sound is obtained by 
the balance of the two equally important (d<j)) A and {dip) 2 terms. There is also a class of 
higher derivative 'galileonic' terms [101111] that can be trusted well beyond the naive EFT 
expectations, as they lead to evolution equations that are not higher order in time. These 
cases are naturally included in our approach. 

• Standard form of the action. In the language of <fi one can always make a field redefinition 
cf) — > <f)(4>). The resulting action is different but it describes the same physics. This redun- 
dancy may be easy to capture for a simple enough Lagrangian but it becomes more difficult 
for a more general one such as, e.g., C(<f), {o~<f)) 2 , □</>)• The action (J4]) is in "standard" form, 
in the sense that there is no residual ambiguity due to possible field redefinitions. 



2 A very general DE theory 

In this section, after a brief presentation of the unitary gauge formalism, we introduce our 
DE action. First we discuss the operators that can be fixed by the background equations of 
motions and then we discuss higher-order terms. 

2.1 Quick guide to unitary gauge 

While referring the reader to [T7l[T8] for more details, it is useful here to give an idea of how 
the unitary gauge works. In a general (perturbed) FRW universe, 4>(t,x) = <fio(t) + 5(j)(t,x). 
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By choosing the coordinate t to be a function of 0, t — t(0), we thus simply have <50 = 0. 
Therefore, as mentioned, the DE action written in this gauge only displays metric degrees of 
freedom. For instance, in eq. ([2]) we have already seen how the kinetic scalar term (<90) 2 is 
transformed in unitary gauge. 

For concreteness, we have made touch with the well known example of a scalar field. But 
what really does for us here is to define a preferred (0 = const.) time slicing. One can well 
forget about and build a Lagrangian with the unit vector n M perpendicular to such a slicing. 
In unitary gauge, 

n = W — S- (5) 

It follows that in the DE action, beside genuinely 4-d covariant terms such as the Ricci scalar 
R, also contractions of tensors with n M are allowed, i.e., by (JSJ), tensors with free upper 
indices such as g 00 , R 00 , etc. Moreover, because time-translations are broken, the coefficients 
of the operators in our action are allowed to be time dependent. These two remarks are both 
exemplified by (J2J). Covariant derivatives of can also be used to express operators in the DE 
action. Equivalently, we can use their projection along and orthogonal to t — const surfaces, 
i.e. the extrinsic curvature 

= h;V a n u , (6) 

where h^ u = g^ u + n^n u is the induced spatial metric, and n a V a n u oc hj^d^g 00 . 

Diffeomorphism invariance is restored by the "Stueckelberg trick", i.e. by performing an 
infinitesimal time diffeomorphism t — > t — ir(x), x — > x. After a trivial coordinate redefinition 
one finds, as a new field in the action, tt(x). This is the Goldstone boson of broken time 
translations, that makes apparent the presence of a new scalar degree of freedom in the DE 
sector. From the action written in unitary gauge, terms containing it are generated by the 
Stueckelberg trick in two ways. First, any time-dependent coefficient, such as c{t) in (j2J), 
generates terms in which tt appears not derivated: 

c{t) c{t + tt) = c{t) + c(t) 71 + -c{t) tt 2 + . . . . (7) 

Second, operators that are not 4-d diffeomorphism invariant transform under time diffeomor- 
phisms and thus generate terms with derivatives acting on tt. Few relevant examples are the 
following, 

^oo ^ g 0O + 2g°> i d f sr + g> a 'd lt ird l >ir, (8) 

SKij ->■ 8K i j + 8 ij d Q {a 2 H)'K-d i dj'K , (10) 
5K -> 5K + 3Htt - a^ 2 V 2 7r . (11) 



2.2 Structure of the DE action: background equations 

Once we assume that there is a metric g^ v minimally coupled to matter fields, we can proceed 
and build the most general unitary gauge action Slg^; t] to describe the dark energy-modified 
gravity sector. In doing this, we can follow step by step the construction outlined in [18] , and 
the general rules there given; with two exceptions: 
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1. We now need to allow a general free function of time f(t) in front of the Ricci scalar. 
In fact, we cannot do any field redefinition of the metric tensor if we want to stick with 
the Jordan frame metric that minimally couples to matter. 

2. For the same reason above, in (J3J) we are not allowed perturbative field redefinitions of 
the metric that re-absorb quadratic curvature terms [42 J such as Ai and A2. 

This procedure univocally determines the only three operators that contain up to linear order 
terms in the perturbations. Schematically, 

+ . (12) 

The universality of the above action is also proved independently by going through the Einstein 
frame construction in Sec. 15.11 and then transforming it back to the Jordan frame (Sec. 15. 2p . 
In unitary gauge S only depends on the metric, which means that we are implicitly allowing, 
beside Einstein gravity, up to one more scalar degree of freedom in the DE sector. 

Action (1 12]) closely resembles a Brans-Dicke theory [SUB]; it is straightforward to recognize 
the meaning of our coefficients /, c and A once we take a Brans-Dicke theory and go to 
unitary gauge. Perhaps slightly less intuitive is to see how more "exotic" operators fit into 
our scheme. While devoting Sec. H]to translate several DE models in our language, here it is 
worth considering, just as an example, the case of an operator for which the correspondence 
with (iT2"l) is less obvious. A natural generalization of the Brans-Dicke coupling <pR, is XR, 
where X = (<90) 2 . Because in unitary gauge X = <filg 00 , and g 00 = — 1 + 5g 00 , the operator in 
question is expressed in unitary gauge as 

XR = <j>t[-R+R( \t) + RW(t)g 00 + 5g 00 5R], (13) 

with R(°) the background value of the Ricci scalar. From the above expression, it is immediate 
to recognize the contributions of XR to /, c and A in equation ffl~2]) . The last term in square 
brackets is quadratic in the perturbations and contribute to the coefficient /i 2 in eq. (J1J). 

It is indeed remarkable that, up to linear order in perturbations, arbitrary complicated 
covariant operators only contribute to the three pieces of eq. (fT2"j) . In what follows we show 
how to fix the coefficients c and A above by using the background equations of motions. In 
the next subsection we discuss quadratic and higher-order operators. 

Let us start fixing c and A by using the equations for a general FRW background of 
constant spatial curvature k. In order to do so, we can neglect quadratic and higher-order 

(2) 

terms contained in S DE . By varying with respect to g^ y we obtain 

UG,u - V M V,/ + g^Uf) Ml + (eg 00 + A)g, u - 2c5°<5° = T, u , (14) 

where T^ u is the energy-momentum tensor of matter. We can take, for instance, the trace and 
the 00 component of the above equation and consider that, on the background, 

G 00 = 3H 2 + 3\, R = -G» = 12H 2 + 6H + 6^, □/(*) = -/- 3#/. (15) 



S 



d X\ 



^f(t)R - A(t) - c(t)g 



00 
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By making the perfect-fluid approximation for matter fields, T" = diag(— p m ,Pm,Pm,Pm), we 
finally obtain, 

M * f 1 ~ k + ^ - \~f + f 7) ~ l^ + Pm) , (16) 

A = Af»/ ( if + 3# 2 + 2 ^2+\j + ~YJ^j ~ \(P™ ~ P») • ( 17 ) 

As already mentioned, the advantage of dealing with Jordan frame quantities is a more di- 
rect connection with observations. For instance, the Hubble constant H = 72 ± 8 km/ s Mpc -1 
measured by the Hubble Space Telescope [35] translates into a statement about the Jordan 
metric (Hubble parameters in Einstein and Jordan frames are related by eq. (|113p ). We also 
refer the reader to the discussion in [27] . In Jordan frame, test particles follow geodesies and 
the energy-momentum tensor in the matter sector is conserved in the usual way, 

p m + 3H (p m +p m ) = . (18) 

In contrast, in presence of a coupling to gravity there is no conserved DE density pn and 
pressure p£>- However, it proves useful to de/in^l the latter through the following (modified) 
Friedmann equations 

H2+ ^ = 3m ipm+pDh (21) 



k 1 

H~ = -TfJ^pipm + PD + Pm + PD) ■ (22) 



Indeed, by using 014p and (|18p . we obtain the modified conservation equation 

p D + 3H(p D + p D ) = 3M 2 J (h 2 + . (23) 



a 2 , 

Finally, using (j2T|)-(|22jh we find 

c = l -{-f + Hf)Ml+ l -{p D+PD ), (24) 
A = !(/ + hHf)Ml + \{p D - p D ) . (25) 

The above are used to obtain the first two lines of eq. (jl]). The only freedom left in fixing c 
and A is thus the function f(t). In practice, however, what discriminates among the different 



3 Equivalently, one could define and p°p as in OH] , such that 

PD = fPD +(f -l)Pm, PD = fpf + (/ - l) Pm ■ (19) 

In this case eqs. (|!?Tj) - (pi2j) get closer to the standard form, 

H ' 2 + i = sip {Pm + pC ° S) ' A - i = ~2^ {pm + p ° + Pm + p ° ] ■ (20) 
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DE models is the choice of few parameters of order one. For instance, by expanding in the 
proper time around the present epoch, 

f(t) = 1 + fMt - t ) + f 2 H 2 (t - t ) 2 + ... . (26) 



In turn, the parameters fi can be constrained by post-Newtonian tests of gravity (see Sec. 
Clearly, an analogous expansion should apply to all time dependent coefficients appearing in 
our action. 

2.3 Quadratic and higher-order operators 

(2) 

The non-universal part of the DE action, S DE , is at least quadratic in the perturbations: 

Sde = f d 4 xy/=g~ F {2 \6g 00 , 6K, U , 6R, upa ; hfV„ t) , (27) 

where is a polynomial of at least quadratic order in the arguments. The first argument 
is Sg 00 = g 00 + 1; hence, includes terms such as 

F (2) D ^2 ( ^00 )2+ M ( ^00 )3+ ; (2g) 

where M 2 4 , M|, etc., are functions of time. In this formalism, a P(0, X) theory looks like an 
expansion in 5g 00 as above (see Sec. 14.11) . 

By considering derivatives projected on the surfaces of constant t, F^ will include also 
terms such as 

F (2) D rot^^oo^oo + _ ( 29) 

Although such terms have not been particularly discussed in the original EFT of inflation 
papers, they are of relevance, for instance, for Lorentz violating models of DE (see Sec. 14.41) . 

The second argument of F^ 2 \ SK^ = — K$ (K^J = Hh^), is the perturbation 
of the extrinsic curvature of surfaces of constant 0, defined in eq. ([6]). We show in App. lAl 
how such tensorial perturbations are defined here. With the addition of this operator and its 
contractions we generate also terms such as 

F m D -^8g 00 5K - ^5K 2 - *3-8K»5K>* u + ... . (30) 

The operators M 2 and Mf are those appearing in the ghost condensate [TH]. We will discuss 
ml, which appear in 'galileon' and braiding models, in Sec. 14.21 

It is worth asking why a term that starts linear in the perturbations such as K does not 
appear in the universal part of fTT2l) . Following appendix B of Ref. [18] we can indeed rewrite 
it as 

(31) 

This relation will be often used in the following of the paper. Note that there is a sign difference 
with respect to eq. (80) of Ref. [18], due to the different convention for the sign of n M defined 
in eq. ()5]) above. 
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Finally, also 4-d curvature invariants are allowed. These are gauge invariant terms and, 
as such, do not produce any n through the Stueckelberg trick. In the absence of matter, it 
is always possible to re-absorb the terms quadratic in the curvature with a perturbative field 
redefinition of g^ v (see e.g. [12] )• However, in the presence of matter this would generate scalar 
couplings to matter fields, which would make us exit the Jordan frame. In App. |A] we show 
how to express curvature invariants in terms of the zero and first-order terms and in terms of 
quadratic operators such as 

d \ x 5R 2 + \ 2 SR^SR^ + lx C^C^ pa + i^ vp °C^ v kX C^ kX + . . . , (32) 

where instead of the perturbation of the Riemann tensor 5Rp UptJ we are using the Weyl tensor 
C^pa [36], defined by 

1 R 

Cpvpa — Rpvpa ~~^S) [l.pRvo QpcrRvp dupRpcT 9vaRpp) ~~^r(.9pp9is<r 9vp9pcr) • (33) 

This vanishes on a FRW background and thus do not need to be perturbed. Therefore, all 
these terms are also included in our action. 



2.4 Comparison with the covariant approach 

As already observed, one of the advantages of our EFT is that only three functions specifying 
the action (jl]), i.e. f(t), A(t) and c(t), are completely fixed by the background evolution. All 
other operators leave the background unchanged and can only be constrained by observations 
involving metric and ir fluctuations, at an increasingly high order in perturbation theory. For 
instance, only the operators in the second and third lines of eq. fl4]) affect the linear evolution 
of perturbations, the other functions entering only at second or higher order. 

This is different from the standard descriptions of DE as a scalar field <fi in terms of a 
covariant Lagrangian. For instance in fc-essence one usually deals with the whole function 
P(4>, (d<f)) 2 ). Although it is elegant to have a covariant description potentially encompassing a 
variety of solutions, an EFT approach should just content itself with the less ambitious task of 
providing a theoretical framework which correctly reproduces cosmological observations, given 
a few parameters. 

EFT can also be directly used at the level of a fully covariant theory. In cosmology, such an 
approach has been advocated and developed by Weinberg in the context of inflation j3H] and 
soon after extended to DE by Park, Zurek and Watson [10] and Bloomfield and Flanagan [12] , 
in Jordan and Einstein frames respectively. The idea is that most models of dark energy and 
modified gravity can be described, schematically, by a scalar-tensor theory with higher-order 
corrections. Schematically, the Jordan frame action is 



S = d A Xy/^g 



M 2 1 

-ffi$)R - -(^) 2 - v(<j>) + m <r\ 



(34) 



What characterizes different models is the specific choice of f(4>), V((f>) and ^[(p, g^], a local 
scalar function of 0, g^ and their derivatives. 

Here we show that, given our unitary gauge action (J3J) with a finite order of operators, 
one can construct a covariant EFT action which correctly reproduces the unperturbed and 
perturbed evolution at that given order (see also [9j[18]). Being covariant, this description 
unifies the evolution of the background with that of perturbations. 
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Let us consider eq. (jlj) up to second order in our expansion, i.e. the first, second and third 
lines of this equation and no other higher-order operator. To make this action covariant we 
perform the following replacements: 




Mt 



sk; - Hit) h; - ^f^<t> 



where R^ and i?/^ are the background values of R and R^ u , respectively given in eqs. ( I119P 
and fll20p of App. |X] The terms in eq. 01]) containing derivatives of g 00 only contribute to order 
higher than fourth in the covariant derivative expansion, and can thus be dropped. Keeping 
only terms up to fourth order in derivatives and integrating by parts several times, it is lengthy 
but straightforward to derive the following covariant action written in Jordan frame, 



By an appropriate choice of functions Oj(^), and Cj(0)Jj this action exactly reproduces 
the background and linear dynamics of (jlj). Moreover, by a number of field-redefinitions it 
can be reduced to the final covariant EFT actions obtained in [T0l[T2] . 

4 The functions ai and bi in the covariant action (|36[) arc related to the parameters of action ^ by 



12Ai(2ff 2 + H) + 2X 2 (SH 2 + H) - m 
~ M 2 
-2V = - 2A + il/ 2 4 + SHfhl - 3ff 2 (3M 2 2 + M 3 2 ) + 36Ai(2H 2 + H) 2 
+ 12A 2 (3# 4 + H 2 + 3H 2 H) - 6/! 2 (2i? 2 + H) , 

c 2M[ + 3Hff4 -6^\(2H 2 + H) mf + 9(#M 2 2 )' + 3(iJAf 3 2 )' 



where a prime denotes the derivative with respect to <f> and all the quantities on the right hand side are now 
functions of 4> and a flat FRW background is assumed. 



s 




(36) 
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3 Scalar dynamics and mixing with gravity 



As already noted, by writing our most general DE action in unitary gauge, we are allowing the 
presence of an additional scalar degree of freedom in the DE sector. In most cases, this can 
be explicitly seen by applying the Stueckelberg trick, i.e., by forcing a time diffeomorphism 
t — > t — 7r(x) upon the unitary gauge action (TjJ, as outlined in Sec. 12.11 

The simplest way to generate a dynamical ir field is considering in our action fll]) a non- 
vanishing c(t). In this case the Stueckelberg trick generates it with a relativistic kinetic 
Lagrangian tt 2 — (V7r) 2 . Another simple example is adding the operator M|. We will consider 
a linearly perturbed flat FRW metric in Newtonian gauge, focussing on scalar perturbations 
only, in which case 

ds 2 = -(1 + 2<$>)dt 2 + a 2 (t)(l - 2V)S ij dx i dx j . (38) 

One way of studying the mixing is considering the linearized equations, which we derive in 
App. IB. 21 In this section we work directly at the level of the quadratic action. By making use 
of eq. OH]) and of the expression for g 00 , one finds 5g 00 — > 2($ — n) + 4$7r — tt 2 + a~ 2 (V7r) 2 . 
Thus, the Lagrangian of it reads, neglecting the expansion of the Universe, 

_ c 5g°0 + M(fy00)2 kin^ic (c + 2M 4 )? . 2 _ _ ^ + ^ $ 

At high energy the last term can be neglected, it decouples from gravity and propagates with 
a speed of sound c 2 = c/(c + 2M|). This is the so-called decoupling limit [18], which takes 
place at an energy higher than E mix ~ (c + M|)/[(c + M 2 4 ) 1 / 2 M*]. 

While decoupling happens in simple cases such as the one above, in this section we will 
study examples where the decoupling is not at work and it and gravity are mixed already 
at the kinetic level. In particular, we will study the mixing produced by i) a non-vanishing 
/; ii) the operator fh\. These are operators that produce mixing but do not involve higher 
derivatives. We will not consider quadratic operators involving higher spatial derivatives such 
as M.2 and M3, which in their decoupling limit have been extensively studied in the context 
of the ghost condensate in [T7HT9] . nor m 2 , studied in [68TT70] . 

Before doing that, let us first mention a well-known case where n is not produced via the 
Stueckelberg trick and yet the theory does contain a propagating degree of freedom. 



3.1 The "hidden" scalar 

By looking at ffl2"]) . it is clear that a way not to generate any 7r is having / and A constant and 
c = S DE = 0, which is the case of ACDM. It is useful to briefly review (following e.g. [9|IT9]) 
why pure Einstein gravity does not contain any scalar propagating degree of freedom. Let 
us expand the Einstein-Hilbert action (see App. IB.lj) and retain only the kinetic terms; more 



specifically, we drop those terms that are subdominant on scales much smaller than Hubble. 
For R, this leaves us with only the terms displayed in the first line of eq. (I134p and we obtain 
(in this section V is a 3-d gradient and V 2 the Laplacian), 

l -V^9R ki = " -3^ 2 + (W) 2 - 2V$W . (40) 

By going to Fourier space and diagonalizing the kinetic matrix corresponding to the above 
expression, we simply obtain the dispersion relation k A = 0, which means that $ and ^ are 
non-propagating degrees of freedom. 
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However, the fact that no n is produced by the Stueckelberg trick does not guarantee, by 
itself, the absence of scalar degrees of freedom. One straightforward example is F(R) gravity. 
By Taylor expanding around a background value R^ its Lagrangian reads 



\y/=gF{R) = l -V=g{F + F'5R + l -F"5R 2 + . . .) 
1 



^(F'R + F- F'RM + -F"5R 2 
2 v yv 2 



+ •••), 



(41) 



where F', F", etc., are computed at R — R(°\ The above equation can be reproduced by 
our action QE) with / = F', A = M 2 (-F + F'RW), c = 0, Ai = M 2 F"/2, etc. Since R is 
a curvature invariant, it does not produce any ir field under a time diffeomorphism, nor does 
any function F(R) thereof. 

However, by rewriting the first line of eq. 04ip at quadratic order in the perturbations using 
eq. (I133p in App. IB. II and taking the high-energy limit, we find 



^V=gF(R) kin ^ tic F' 



(V*) 2 -2V$W 



3^ - 2V 2 * + V 2 $ 



n 2 



Calculating the determinant of the Lagrangian in Fourier space one finds 

det £ = k 4 F' [3F"(u 2 - k 2 ) - F'] . 



(42) 



(43) 



Setting this determinant to zero gives the dispersion relation of a relativistic degree of freedom 
with mass m 2 = F'/3F". 

More explicitly, we can consider the constraint equation 5S/5<& = derived from eq. fj4*2|) . 
Every term of this equation is acted upon by at least one Laplacian V 2 that can be safely 
integrated, giving 

F' 



F" 



= 2V^ 

Plugging this back into the action, integrating by parts we finally obtain 



-g_ F{R) q uad 3F/ 



- (v^y 



F" 



(44) 



(45) 



the action of a relativistic scalar field. The condition to get to (]45p is that F" ^ 0. In 
Sec. 14.51 we outline an equivalent way to deal with F(R) theories which is closer in spirit to 
our formalism and brings the action to the universal form (I12p with Sp E = 0. 



3.2 Mixing proportional to / 

We now turn to the situation in which a tt field is produced via Stueckelberg trick and it is 
kinetically mixed with gravity. Here we consider the well-known case of the mixing produced 
by a non-constant f(t). This is the type of mixing common to all scalar-tensor theories and 
we just review it here briefly in our language. By using the results of App. \B.1\ we write the 
relevant kinetic terms in Newtonian gauge, 



S -i 
b ~ 2 

kinetic 



M 2 f [-3* 2 - 2V$V* + (W) 2 + 3(///)*7T - (///)vr(V 2 <l> - 2V 2 ^) 



(46) 
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If we go to Fourier space and calculate the determinant of this 3x3 kinetic matrix we find 
~ k 4 (u 2 — k 2 ), signaling the presence of a single propagating degree of freedom with relativistic 
dispersion relation uj 2 = k 2 . 

The transformation that de-mixes tt from gravity is well known, it is the conformal trans- 
formation to the Einstein frame metric, see Sec. [5j On the metric potentials $ and \& this 
transformation reads 



and action fHBTl becomes 



M 2 J 



-3*1 



2V$ B V* B + (W E ) 2 + -(///) 



2 

71 



(Vtt) 



(47) 



(4? 



Even though in this frame it is kinetically decoupled from gravity, it is now directly coupled 
to matter. In order to study the effects of this theory on matter we consider its Newtonian 
limit: we can neglect time derivates and introduce pressureless matter sources!! We will do 
that directly in Jordan frame because this is the frame where test particles follow geodesies. 
For generality, we also add to eq. (l4"B"j) an explicit tt kinetic term — c(dir) 2 . Varying this action 
with respect to * gives 

* = $ + (///> • (49) 
Plugging this back into (I46I) finally yields 



5 



M 2 J 



-(V*)' -(///) V* Vtt - (///) 



M 2 f 



(Vtt) 



<$>6p r 



(50) 



where we have added the gravitational coupling to non-relativistic matter (l/2)T 00 5g 00 . 
Varying this action with respect to tt fixes 



TT 



-Mlf 



which can be used in eq. 
and floTj) we get 



2(c + Mlf 2 /f) 
to derive the post-Newtonian parameter 7 = */$. From 



(51) 



1-7 



Mlf / f 



(52) 



2(c + Mlf 2 /f) 

By varying with respect to $ we see that the usual Poisson equation for the gravitational 
potential is modified by a term in tt. Finally, using (15 ip the Poisson equation reads 



Mlf / f 



4(c+M*P/f) 



V 2 $ 



1 



2Mlf 



5p r 



(53) 



By comparison with V 2 $ = p m /2Mp l we finally obtain the relation between M* and the 
measured Newton constant in the presence of a fifth force 



1 c + Mlf 2 /f 
MPJ c+\MlP/f- 



(54) 



5 We discuss the Newtonian limit in the quadratic approximation. This means that we are neglecting the 
possible screening effects due to the higher-order operators. 
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3.3 Mixing proportional to ml 

Here we consider in Sp E of eq. (I12p the operator which can also produce kinetic mixing [9]. 
This operator is typically contained in galileon [^QlSllSTlSH] and braiding models [19], as we 
discussed in Sees. 14.21 and 14.31 We can read off the effect of the Stueckelberg trick on this 
operator from eqs. (jHJ) and (TTTT) . In particular, in Newtonian gauge at linear order we have 



5g 00 ->■ -2(tt - $) , c^-> -(3* + a" 2 V 2 7r) . 



(55) 



Thus, by keeping only the kinetic terms the action reads 
S = 



kinetic 



M. 



r2 



2 " 2 

-3^ 2 - 2V$W + (W) 2 



(56) 



C7T 2 



c(Vvr) 2 - 3m^7T - mfV^Vn 



Here, as before, we have dropped the scale factor since we are considering scales much smaller 
than Hubble. Also, we have defined 



c = c+ -(Hm\ + m\) 



(57) 



where the second term on the right-hand side of this equation comes from integrating by parts 
the piece — amf7rV 2 7r. 

It is handy to re-define the fields in such a way that they all have dimensions of mass as 
follows, 



7T C = C 1/2 7T , ^ c = , $ c = M*$ 

The Lagrangian becomes 



vrci 



a 



2M,cV2 



£ = -3^ 2 - 2V$ C W C + (W c ) 2 + 7T 2 - (g/c)(V7r c ) 2 - 6a^ c 7T c - 2aV$ c V7r c 
In Fourier space, on the (\l/ c $ c 7r c ) basis, the Lagrangian is a 3 x 3 matrix, 

-3co 2 + k 2 -k 2 -3au 2 \ /*A 



(58 



(59) 



£ ~ (tf c $ c 7T C ) | -fc 2 

-3aw 2 

The determinant amounts to 



-ak 2 
-ak 2 io 2 -(c/c)k 2 J 



(60) 



det£ = -k A [u 2 (l + 3a 2 ) - k 2 (c/c - a 2 )] 



(61) 



which corresponds to having one propagating degree of freedom with dispersion relation u 2 
c^k 2 , where the sound speed c 2 is 



m\/Ml 



c + \m\jMl 



(62) 



Note that c 2 can become negative, in which case one has a gradient instability. We will get 
back to this point when we study the stability in Sec. 13.41 
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One may wonder whether there exists a transformation of the metric that de-mixes it from 
gravity, as in the example of the previous subsection. Indeed, it is easy to show that the two 
new potentials 



1 ™- lTr lTr I J- 



= $ + — i-vr, ^ E = ^ + — i-vr, (63) 

kinetically de-mix the system. Comparing with eq. (I47p . one sees that this cannot come from 
a conformal transformation of the metric. 

As before, we can now add matter and study the Newtonian limit. We can neglect time 
derivates and, by varying the action (1561) with respect to obtain $ = That the two 
potentials are equal at linear order in the absence of anisotropic stresses is also evident from 
the traceless ij part of the linearized Einstein equations, eq. (j!49p of App. IB. 21 This is an 
example of a modification of gravity where, even in the absence of screening and in the presence 
of an explicit kinetic mixing, the post- Newtonian 7 parameter is unity. In this case, the action 
fl56l) becomes 



S = J ' -Ml(V§) 2 -£(Vvr) 2 -m^WVvr - <Mp m . (64) 

Analogously to what we did in Sec. 13.21 we derive the observed Newton constant, 

r = J c+\(rn\ + Hm\) 

N 8ttM2 c - \m\jMl + \ (rn\ + Hm\) ' 1 ' 

where we have used the definition of c, eq. (15T|) . We stress again that the above analysis holds 
in the presence of a long-range propagating scalar degree of freedom and we have not assumed 
any screening mechanism. 

3.4 De-mixing in a quite general case and stability 

Now we turn to a more general situation in which several of the discussed operators are 
switched on: 

S = j (j£fR - A - eg 00 + ^(6g°r - ^-5g 00 6K^ + \^gT^5g, v . (66) 

Here, in the last term, we have also added the coupling of matter to gravity. This will give 
us the opportunity to discuss the effect of fifth- force mediation on the Newton constant. Note 
that for the present discussion A is irrelevant. 

Once a DE model is translated in our framework and "distilled" into a certain number 
of operators, it is straightforward to address its stability by de-mixing it from gravity and 
looking at its quadratic action. Only / and fn\ produce kinetic mixing in the above action. 
We should then combine the linear transformations fT4T|) and fl63|) and define 

$ R = $ + - \ J - + — h 7T 



2\f Ml 
Iff fn\ 



(67) 
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The kinetic quadratic piece of action fl66|) written on this new basis features a normal non- 
propagating gravity sector (see eq. (HOI) ) plus the following de-mixed scalar piece: 



kinetic 



2M 2 4 



3p 
4/ 



3- 3 / 
-rrti — 
2 V 



3 

4M? 



7T 



1 3 / 2 , ,2 1 

c H M," ,i>, — 

4/ * 2 V 



3/ _ 1^1 + 1 ^3 + 7^3\ 

4 AT 2 2 V 1 



(Vvr) 



• (68) 



The condition that there are no ghost instabilities requires that the brackets in the first line of 
this equation is positive. Note that by eq. ( 12~4"|) c = (pn + pu)/2 — (f — Hf)M%. Thus, by an 
appropriate choice of parameters the equation of state may become phantom, po < ~Pd, and 
yet be ghost free — for instance, when M| dominates the kinetic term [9j[T7], in the presence 
of a time dependent / [50J or of a non- vanishing m\ jl9j . 

The square of the speed of sound of fluctuations is given, as usual, by the ratio of the 
brackets in the second line of action (1081 to the one in the first line, i.e., 



c +lMlP/f-\m\f/f 



1 m\/Ml + \{fn\ + Hm\) 



c + 2M| + \MIPU - \m\f/f + \m\jMl 



(69) 



Since the denominator on the right-hand side is positive by the ghost-free condition, c 2 has 
the same sign as the numerator. This may become negative, in which case one has a gradient 
instability. In particular, regardless of the sign of m\ the mixing with gravity always induces 
a Jeans-like instability proportional to m\. This instability is discussed in [9j[T7], where it is 
assumed that M| dominates the time kinetic term and that / and m\ are constant, in which 
case it is cured by requiring that m\/Ml < H . More generally, ensuring that the time for these 
instabilities to develop must be longer than the age of the Universe constrains the parameters 
of action (HJ. 

The calculation of the modified Newton constant for this DE theory can be done straight- 
forwardly using the methods developed in the previous two subsections. It yields 



G, 



c + MtP/f + \(ml + Hm\) 



SnMlf c + *M*P/f - \m\f/f - \m\jMl + \ {m\ + Hm\) 



(70) 



4 Matching with DE models 

To show that the action (j4j) is very general and encompasses many of the models of the 
literature, here we discuss how one can translate in the EFT language several well-known 
examples. 



4.1 Quintessence and /c-essence 

For f{t) = 1, the action (j3J) contains any minimally-coupled single-field dark energy model [0]. 
For instance, as discussed in the introduction, a model with minimal kinetic term such as 
standard quintessence is contained in its first line. Indeed, its Lagrangian can be rewritten in 
unitary gauge as 

- kd<t>) 2 - V{4>) ~Ulit)9 m ~ V(t) , (71) 
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which is reproduced by the first line of (jlj), once we note that 2c(t) = po + Pd = <Po(t) and 
2A(t) = p D -p D = 2V(t). 

Similarly, the action (jlj) includes also /c-essence models jST], where the Lagrangian has at 
most one derivative acting on the field, £ = P(<f),X). In unitary gauge this is of the form 
P(4>o(t), 4>lg 00 ), which can be expanded in powers of ipo^Q 00 an d written in the form (jlj) [9|H~8]. 
Note that one can always make a field redefinition such that 0o — t. In this case the action 
(jlj) is reproduced with 



A(t)=P(t), c(t) 



dP 
dX 



M 4 n (t) 



gap 



X=l 



[n > 2) 



(72) 



x=x 



4.2 DGP and galileon 

Terms containing the extrinsic curvature can become important around localized gravitational 
sources and be responsible for screening mechanisms, such as the Vainshtein effect [61] , which 
allow models based on light fields with gravitational-strength coupling to evade Solar System 
constraints. In particular, in the decoupling limit the DGP model [65] is described by a 
non- minimally coupled scalar field with cubic self- interaction [17] H8]. 



m° 



(73) 



which is the prototypical operator responsible for the Vainshtein mechanism and which is also 
the cubic operator £3 of the galileon [10]. Now we want to rewrite this operator in unitary 
gauge to show that it is actually contained in action (jlj). Specifically, we can expand the 
operator (I75|) in powers of 5g 00 and at most one power of 5K (see also [66]). 

Using the definition of the extrinsic curvature, eq. (jBj), we note that, up to a boundary 
term, the density Lagrangian can be rewritten as 



(d<pfn<p=^-[-(d<pf} i/2 K . 



(74) 



In unitary gauge this becomes 
(90) 2 D0^(2/3)(V°) 3/2 0^ 



0o - H<p ((V°) 3/2 - 1) + (2/3) [(- 5 00 ) 3/2 - 1] 



(75) 



where for the second line we have used K 
(2/3)0 3 if using 

K = -V„ 



3H + 5K and we have integrated by parts 



//O 



no 



(76) 



Thus, the DGP operator (173j) can be written in the EFT language as a polynomial of Sg 00 and 
SK with at most one power of SK. 

As a concrete example of a model which contains the operator (1731) . let us consider the 
galilean cosmology discussed by Chow and Khoury in Ref. [67] and rewrite it in our EFT 
language. The full action in Jordan frame reads (see eq. (3.1) of this reference) 



S 



d A Xy/—g 



M 2 r 2 
^± e -24>/M tR _ ^ W) 2 D0 + c 

2 M* 



(77) 
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where r c is a length-dimension parameter. We can expand this action in series of 5g 00 and 5K 
by making use of eq. (T75|) . Comparing with (J3J) we find 

/(t)=e- 2 ^ , A(t) = _^02(0 o + 3#0 O ), c(t) = -^02 (0o _ 3 ^ 0o)! (7g) 
M^(t) = 2 (0o + 3#0 O ) , M 3 4 (t) = -^0o(0o + #0o) , M 4 4 (t) = . . . , (79) 

2 

m?(t) = -^-20* , m^) = ..., M 2 (t) = M 3 (t) = ... = 0. (80) 

As an application of the machinery of Sec. (13.41) . it is possible to study the stability of this 
model. Our analysis is slightly different from that of [67], as it takes into account the kinetic 
mixing with gravity. Using eq. (1681) . the requirement that there is no ghost instability turns 
into 

~ 2 M 4 "' H + W^- 2 W^ + e ~^ >0 ' < 81) 

while the speed of sound squared, replacing eqs. fl78"l) - fl8"Ul) into eq. (I6"9~|) . reads 

c2 = -IHM*/M* - \igrHMl - \4>yjMi ~ IM 2 JM* + 
-2#0 o r c 2 /M* - 2$rl/M* + 4 r c 4 /M, 4 + e -^°/ M * 

Using the approximate solution studied in O ~ ±(2/3) 1 ' 2 M*r~ 1 , H « (2/3) 3//2 r~ 1 and 
0o *C M*, we confirm that the stable branch is for O < 0, for which the sound speed squared 
( 182]) exactly reduces to c 2 s = 1, a different result from that of [67], i.e. c 2 s = 2/3. 

4.3 Kinetic braiding 

It is straightforward to show that the operator 

G(0,X)D0, (83) 

denoted in [19] as kinetic braiding, can also be rewritten as a polynomial of Sg 00 and SK with 
at most one power of 5K. We can expand G(4>,X) in polynomials of SX, say, 

oo 

G(0, X) = V^X U<P)$X m , (84) 

m=0 



where the y—X on the right-hand side helps simplifying the notation and Z m (0) are the 
coefficients of the expansion. Then we can use the relation 

1 



□0 = -V^XK - ——n^X . (85) 
2v — X 

In unitary gauge X = 0q(— 1 + Sg 00 ) and we will take 0o = 1 for simplicity. Terms linear in K 
can be integrated by parts as in (|3"T1) and give powers of 5g 00 . After few integrations by parts 
we finally obtain 



-OOym 



G(0,X)D0 = k(i) V C^ + £(_ /m (t)ff + ?[^ 

(86) 



m=l 
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The contributions of such an operator to action (j3J) read 



A(t) = 3H^h-^l J , c{t) = -i + 3H^h-^l J (87) 

M 4 2 {t) = -j + | - 61H 2 + ^fT/i, M 3 4 (t) = . . . , (88) 
m\{t) = 2h - 3/ , rhl(t) = ... . (89) 



4.4 Ghost condensate and khronon 

Operators containing 5K 2 , (5K V ^) 2 and gradients of g 00 have more than one derivative acting 
on a scalar. These terms can become relevant in models where the spatial gradients are 
parametrically more important than the time derivatives. 

For instance, the ghost condensate theory [19], based on the field shift symmetry — > 
0+ const, is realized in the limit of (pd+Pd) -+ 0, i.e. c(t) — > — here we are assuming minimal 
coupling, / = 1. In this limit the spatial kinetic term of the perturbations is dominated by the 
higher-derivative operators proportional to M 2 and M 3 [9^ ITT1 - [T9"] . Controlled models with very 
small speed of sound [MHM] should be thought of as small deformation of the ghost condensate 
limit. In this context, absence of ghost and gradient instabilities can be guaranteed also for 
an equation of state po < ~Pd by the presence of these higher-derivative operators [17], even 
though their effect is practically absent on cosmologically relevant scales [9]. 

Another example is the khronon field, based on the full reparametrization invariance (f) — > 
</>(</>). This symmetry has recently received attention in the context of Horava gravity and its 
healthy extensions [68TI70] and it has been used to construct a technically natural dark energy 
model in [7T] and test Lorentz invariance of dark matter [72]. The invariant object under this 
symmetry is the unit vector perpendicular to the constant time hypersurfaces n M — defined in 
()3]) — and to lowest order in derivatives the action can be written as [7TI|[73] 

S = l - J d A xV=g (MlR -2A- M A 2 (V/ - m) 2 + M^KV,/) 2 ) , (90) 

where M\ and M a are the two mass parameters of the model, beside the vacuum energy A. 
It is straightforward to show that in unitary gauge this action is equivalent to (j4j) with_| 

A(t) = A , c(t) = , 

M„ (91) 
M 2 (t) = M A , m2 (t) = -f. 

Thus, symmetry under reparametrization <fi — > <fr(4>) imposes that all the terms without deriva- 
tives on g 00 are absent, i.e. c(t) = Mi(t) = 0. 



4.5 F(R) gravity 

Curvature invariants offer an interesting application of our formalism. We should distinguish 
two types of situations. At the low energy effective level, the gravitational Lagrangian does 
include a series of terms made of curvature invariants, such as R 2 , R^R^ etc., weighted by 

6 To rewrite the last term, use eq. (73) of Rcf. j!8j . 
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the smallest masses that have been integrated out (see e.g. the nice review [12] )• These terms 
do not generate new degrees of freedom and their natural size is generally way too small to be 
relevant for the DE problem. Our formalism can encode curvature invariants; as an example, 
in App. IA1 we consider the covariant operator A{(j))R^ v R IJ,u and calculate its contributions to 
the various terms of action (J4]). 

When higher-order curvature invariants are used, naively, beyond the EFT regime, they 
normally lead to ghost instabilities [52]. As far as we understand, the only exceptions are 
general function of the Ricci scalar R [55H55] and of the Gauss-Bonnet term [55]. Leaving 
aside naturalness problems, there has been a lot of activity in trying to engineer DE models 
using Lagrangians of the F(R) type (see e.g. [57] for a review and further references). The 
equivalence between F(R) and Brans-Dicke theories is well known (see for instance [58]). Due 
to the strong coupling of the resulting scalar-tensor theory, the Chameleon mechanism [59] is 
needed to make F(R) models phenomeno logically acceptable. 

Let us see how to treat F(R) theories within our formalism. First we should note that a 
generic expansion in SR is allowed by the set of rules that we explained in Sec. 12.31 Indeed, 
by expanding around a FRW solution with Ricci scalar R^(t), we obtain the second line of 
eq. (|4"TI) . This is already a legitimate DE action in unitary gauge: as discussed in Sec. (13. ip . 
at quadratic order it generates a propagating degree of freedom with relativistic dispersion 
relation. However, if we attempt the Stueckelberg trick at this point, we find that no tt is 
indeed produced. This is not surprising: neither R, nor any other curvature invariant produce 
7r field upon the infinitesimal time diffeomorphism t — > t — ir(x), just because they are gauge 
invariant. 

However, we can fix the time slicing in such a way that it coincides with uniform R 
hypersurfaces. Note that this does not mean that we are killing the dynamics, because 5R 
has a nontrivial structure when written in terms of metric components. For instance, the 
linear term in SR in the first line of eq. (l4~Ti) still produces non-trivial equations of motions 
when varied with respect to some metric component. What this prescription does kill are all 
the terms beyond the linear order in the expansion (1411) . because their contributions to the 
equations of motion has always at least one power of 5R, and therefore vanish. Finally, by 
choosing the time such that R^ = t from eq. (j4ip we obtain the unitary-gauge action 

/M 2 f M 2 

d A x^-g^-F{R) -> J d'x^-g^- [F'(t)R + F(t) - tF'(t)} . (92) 

This action is equivalent to (jlj) upon identifying 

f(t) = F'(t) , A(t) = ^ (F(t) - tF'(t)) , c(t) = , (g3) 
Mi = Mi = . . . = . 

Since now the gauge is fixed, it terms are produced by expanding the time dependent functions 
F(t) and F'(t). At first sight, tt is non-dynamical because it does not have a kinetic term. 
However, the latter can be generated by de-mixing the field from gravity, as explained in 
Sec. 13.21 Indeed, as expected, action ( 192]) is the one of a scalar field without kinetic term in 
Jordan frame, written in unitary gauge and the above procedure is alternative and equivalent 
to the standard trick that introduces a scalar field <fi in the action [58] . 
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5 The "Einstein Frame" 

The constructive approach explained in the introduction naturally leads to a Jordan frame 
formulation of DE in the unitary gauge. The Jordan metric has also the advantage of being 
uniquely defined by the coupling to matter and more directly related to observables as we noted 
in Sec. 12.21 When / is not constant, however, the gravitational dynamics looks obscure in the 
Jordan frame, just because the new scalar degree of freedom, it, mixes with the gravitational 
ones — say, $ and ^ in Newtonian gauge. In the absence of kinetically mixing quadratic 
operators such as mf , the system can be diagonalized in a very standard way by a conformal 
redefinition of the metric tensor, i.e., by going to the so called Einstein frame, which is the 
subject of this section. However, one should keep in mind that, in the presence of mixing 
quadratic operators, gravity cannot be de-mixed in this way, although it is still possible to 
diagonalize the quadratic kinetic action as shown in Sec. El 

In this section we adopt the standard notation of scalar-tensor theories and call "Einstein 
frame" the conformally related metric in which the gravitational action does not contain a 
function multiplying the Einstein-Hilbert term and in which there is no mixing proportional 
to /. Instead of just conformally transforming all our operators and coefficients, in the next 
subsection we provide an alternative, "Einstein frame" construction of our formalism as a 
natural extension of the EFT of quintessence developed in [9]. We then re-switch back to 
Jordan frame in Sec. 15.21 

5.1 EFT of DE in Einstein frame: alternative const ruction@ 



It was shown in [18] that the most generic theory with broken time diffeomorphisms around 
a given FRW background can be written as 



S = d A x\/^g 



2 

where c and A are fixed by the background evolution 



_ cg oo _ a + F^ 00 , SK^, 5R, vpa] V*V M ; t) 



(94) 



c(t) = -Mi (If- -J , A(t) = Mi \3H 2 + H + 2-J , (95) 
and 

is a quadratic or higher-order function of its arguments, which generates terms 
invariant under spatial diffeomorphisms and of quadratic or higher order in the fluctuations 
around a given FRW background |17j . 

Action f )94p was used to describe the EFT of inflation in the context of single field models. 
One of the main differences of the current phase of cosmic acceleration with inflation is the 
presence of matter species such as cold dark matter, baryons, radiation and neutrinos. In 
this case, one must supplement eq. (I94|) with an action for these components. Neglecting 
mutual interactions between components, which are irrelevant for late-time cosmology, the 
action reads 



,5' 



/" n A/f 2 1 
<*V=0 - eg 00 - A + . . . + ^1 , (96) 



7 This subsection and App. IC.ll are the only places in the paper where Einstein-frame quantities are not 
denoted with a hat nor with the subscript E and Jordan-frame quantities explicitly carry a subscript J. 
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where the index i runs over all matter species and the ellipses stand for higher-order terms. 
The functions c and A can be fixed imposing that the background energy density and pressure 
of the dark component p£> and pr> satisfy the Friedmann equations. More specifically, the 
background Einstein equations derived from this action read 

G^Ml + (A - c)g, u - 2c$$ = ^ T« . (97) 

i 

Thus, using (1151) in this equation yields the Friedmann equations with 

<t) = \ {Pd + Pd) , A(t) = \{Pd- Pd) ■ (98) 

Once the linear terms have been fixed using ( 1981) . action (19"6T) describes the most generic theory 
of single-field dark energy perturbations around a give FRW background evolution. We stress 
that in this action the higher-order terms in the dots are constructed from the metric and from 
geometrical objects describing the preferred time-slicing and do not explicitly involve matter 
quantities, such as the matter energy-momentum tensor. Thus, matter is minimally coupled 
and for this reason this theory has been dubbed EFT of quintessence [9]. 

To go beyond the assumption of minimal coupling, we can couple each matter species to a 
metric gffi which differs from the Einstein frame metric, now denoted by g^, for which the 
mixing between dark energy and gravity is minimal. Schematically, 

SmAg^M — ► Sm^gW,^] . (99) 

For each matter species, the metric gfiu ^ is a two-index tensor whose most generic form 
can be given as a function of the Einstein frame metric g^ u and of geometrical objects in- 
variant under spatial diffeomorphisms, constructed in the Einstein frame. We will assume for 
simplicity that all species couple to the DE sector in the same way, but our arguments can 
be straightforwardly extended to considering non-universal coupling^]. In unitary gauge the 
Jordan metric takes generically the form 



msg oo +lm(sg°r 



Sg°^g 0u [ 7 o(t) + 7i(W° + ...] 



+ ^\t)g^6g 00 K + Knf{t)8K^8K; + . . . , (100) 

where f(t), (3i(t), ji(t) and rriiit) are functions of time. The first line of this equation includes 
conformal and disformal transformations |74j . allowing for kinetic mixing. Note that the 
Jordan metric has been written in such a way that the only term contributing to its background 
is the first on the right-hand side, f{t)g^ u \ the only linear term is the one proportional to Pi{t), 
while the other terms are higher order in the perturbations. Their exact form is irrelevant 
for our discussion. Indeed, in App. IC.ll we show that all the terms on the right-hand side 



8 The weak equivalence principle is so well constrained (the level is < 1CP 13 [75]) that its possible tiny 
violations are hardly relevant for DE — this is why it is a good approximation to assume a universal coupling 
here. On the other hand, couplings to the dark matter sector are far less constrained (order 10% [751177] ) 
and play a crucial role in models of coupled quintessence (e.g. [31h33] 1 It would be interesting to extend our 
formalism to those scenarios. 
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of fllOOp except the first can be reabsorbed by a perturbative field redefinition, so that, for 
what concerns the matter action, we can writqj 

S m [gTM=S m [f{t)g^^ m \. (101) 

This field redefinition leaves the function A unchanged, while changing c and the time- 
dependent coefficients of the higher-order terms. However, the overall structure of the gravi- 
tational/DE sector action remains the same as eq. In summary, 



S = I d 4 xy/^g 



RM * 00 A 

- co - A 

2 y 

+ ^W°) 2 + fW°) 3 + • • • - ^ Sg m 5K - ^6K 2 + ... 



(102) 



where in the second line we have explicitly written some of the higher-order terms. 

Before deriving the Jordan frame action, let us briefly review the background evolution 
equations in the Einstein frame. In this frame, the background Einstein's equation is not 
affected by the presence of the conformal function f(t) inside the matter action and it is still 
given by 

G^Ml + (A - c)g^ - 2c5j5° = . (103) 

Thus, eq. ( |98|) still holds. However, the background matter energy-momentum tensor in the 
Einstein frame is related to the one in the Jordan frame by T^ u = f~ 2 T^ u , which implies 
Pm = Pmf~ 2 an d Pm = Pmf~ 2 - The Friedmann equations are obtained from (11031) and are the 
usual ones, 

H2 + V 2 = mi {pm+pDh (104) 

k 1 

H ~ ^ = -7^piPm + PD +Pm.+PD) ■ (105) 

As expected, matter and dark sector are coupled in this frame and their conservation equations 
read 

If If 
p m + 3H(p m +p m ) = -jT m , PD + 3H(p D +p D ) = --jT m . (106) 



5.2 Going to Jordan frame 

To go to the Jordan frame we perform a conformal transformation ds 2 E = fds 2 j in the action 
(11021) . In this subsection, contrarily to what we did in the previous one, we drop the index J 
from Jordan frame (or "J-frame" , for brevity) quantities and we use the index E for Einstein 
frame (or "E-frame") quantities. At times, E- frame quantities are also denoted by a hat to 
make the notation more compact. We have 

d A x E ^=g^ = d 4 xf 2 ^ , (107) 

«-r>(«-3nM-5 9WW ) . (1W) 

9 A similar field redefinition to rewrite the action in a minimal form has been also done in [12) . see e.g. their 
equation (1.2). 
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Things simplify by rescaling the time and the scale factor such that (He — f l ^ 2 dt and a# = 
/ 1 / 2 a, while leaving spatial coordinates unchanged. With this choice, in unitary gauge the 
conformal transformation of the metric components readJ^I 



9oo = 9oo, 9oi = (a E /a)g 0i 



9 



ij 



{a E /a) 2 g i: j 



(109) 



and simplifies the transformations in the unitary-gauge quantities in the action. 

It is convenient to first discuss how the first and last lines of the action (11021) transform, 
discarding for the moment the higher-order terms in the second line and leaving the discussion 
of their conformal transformation to the end of this section. Indeed, as expected, transforming 
the higher-order terms does not affect the first and last lines of the action. In this case, the 
action in Jordan frame reads as in eq. (1121) . with c and A related to the analogous E-frame 
quantities by 

3 fd\nf s 

4 I dt E 



f 



ce 



Mi 



A = f 2 A E 



(110) 



The J-frame background energy density and pressure defined in (l2~Tj) - f[2"2l are related to their 
E-frame counterparts introduced in the last subsection by 



PD = f 



PD = f 



^ din/ / l dln/ 

dtE \4 dtE 



Pd + 3M 2 



-Hr 



Pd + Ml 



d 2 In/ 

dt\ 



din/ 

dtp. 



dlnf 



dtj 



112) 



Note that the Hubble parameter is different in the two frames, 



f 



-1/2 



H+ 2j 



;ii3) 



where a dot denotes a derivative with respect to the J-frame time. 

For the higher-order terms in the action, let us start by discussing the polynomials of 
5g°°. Since Sg 00 is not affected by a conformal transformation in unitary gauge (see eq. (I109p ). 
the coefficients of this expansion will be simply multiplied by the factor f 2 coming from the 
transformation of the volume element, eq. (11071) . Thus, in Jordan frame the first two terms of 
the second line of eq. (11021) simply become 



f 2 M 4 f 2 M 4 

J M * (6g 00 ) 2 + L3L( S g 00 ) 3 + 



2! 



3! 



(114) 



We can now discuss the terms containing the extrinsic curvature. From the definition of the 
extrinsic curvature and using eq. (I109p . in unitary gauge one obtains 



K» v = f 



-1/2 



-9 00 f 



f 



(115) 



10 Note that these relations only hold in unitary gauge. Some care has to be taken when extending them to 
an arbitrary gauge. For instance, the conformal transformation ds 2 ^ = fds 2 j in unitary gauge reads g^dt 2 ^, = 
fgoodt 2 , which in a general gauge becomes, at linear order, g^ dt 2 E = /(l + (f/f)ir)goodt 2 . In Newtonian 
gauge this yields $b = $ + (///)tt, which is of course consistent with eq. (|T7)) . 
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so that 



= r 1/2 1 5K» - - S9 °° -h" 



(116) 



where 1 / \/—g 00 can be expanded in polynomials of 5g 00 . Including in the conformal transfor- 
mation the two last terms containing 8K, the second line of eq. f l 1 2 j) becomes 



d 4 xy/—g 



t2 £2 



+ 



2 8 32/ 

PMI , 3/V^ ,,. U |/2 



3! 16 

pl^ml 3M 2 2 / ' 



64/ 




117) 



From the above, one can now easily recognize all the different pieces in the J-frame action (j3J) 
that we re-derived in the last two subsections from the E-frame perspective. The inclusion 
of other terms in the transformation is straightforward as well as the inverse transformation 
from Jordan to Einstein frame. The latter is done explicitly in App. IC.2I 



6 Concluding remarks 

We proposed a description for dark energy and modified gravity that aims to encode all models 
in a unifying fashion. Our DE action (T4|) is organized as an expansion of increasing order in 
number of perturbations. This is made natural by the use of unitary gauge. In this gauge 
the dynamics of both the gravitational and scalar sectors is encoded in the degrees of freedom 
of the metric. Therefore, our DE action is written directly in terms of the perturbations of 
the metric field around a FRW solution. Only three operators, fit), c{t) and A(t), contribute 
to linear order and are thus fixed by the background evolution. All the others start at least 
quadratic in the perturbations. Each operator is responsible for distinctive dynamical features. 

In Sec. [3] we initiated a systematic study of such features. We considered a restricted 
number of notable operators in the limit of high energy. This is sufficient to study the ghost 
and gradient stability and to compute the speed of sound of the associated scalar fluctuation— 
we leave a more detailed analysis of the Jeans instability and of the clustering behavior to 
future work. Note that all those features can be studied once and for all in our formalism, and 
then any new DE model is analyzed straightforwardly, once it is "distilled" into our basic set 
of operators. For example, once the coefficients /, c, M| and rh\ have been worked out for a 
given DE model, its stability and speed of sound will be simply given by eqs. (168)) and ( 1691 ). 

Our formalism also allows to address general issues in an efficient and easy way. For 
instance, there has been interest in understanding whether one can have a sensible theory with 
violations of the null energy condition or, in other words, a "super-accelerating" equation of 
state w < — 1. For minimally coupled models this has been addressed in [SJCT7], and we can 
review the argument here by considering the limit / = of our theory. By the background 
equation fT2"4"j) . w < — 1 implies c < 0. We can ask whether the Lagrangian for the fluctuations 
can still have the "right" signs and the theory be stable. By looking at eq. fl68|) it is immediate 
to realize that ghost stability — related to the sign of the n 2 term — is guaranteed for high enough 
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values of M|, fn\. In turn, gradient stability will set an upper limit on the ratio fh\/{M^H) 
(see, e.g., [IH] for an explicit model displaying this property). The present framework allows 
to generalize this argument to non-minimally coupled models (/ ^ 0), which are also known 
to provide an effective equation of state w < — 1 [50]. More generally, the question of the 
stability can be immediately addressed from eqs. (124)) and (168|) . where reduces to an algebraic 
problem. Without need of tedious model building: a large variety of scenarios can be explored 
simply by varying in different ways the relative contributions of c, Mf / and fh\. 

Our description looks handy also for an efficient comparison against the data. The main 
reasons are the formulation in terms of an action principle and the fact that effects of quadratic 
and higher-order operators are disentangled from the background evolution. A more explicit 
link with observations is left to future work. In particular, it will be important to extend the 
dynamics beyond the strict high-energy limit considered here, in which case matter species 
acquire an important role. The aim is to match the operators of our description to observation 
of the large-scale structures, both in the linear and non-linear regime, and including higher- 
order statistics. 
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A Curvature invariants 



Here we explain how to include in the unitary-gauge action general curvature invariant terms 
by considering R flI/ R' IU as an example. We follow closely (and refine here and there) the 
method outlined in App. B of Ref. EH 



All curvature invariants are made with the Riemann tensor R 



pupa • 



(2") 

In S DE we want to 



include terms that start explicitly quadratic or higher in the perturbations, which means 
polynomials of at least second order in 5R^ upa = R^p^— Rpupa, 5R pu = R pv — Rpu or 5R = 
R-R(°\ The trick used in |18j is to exploit the high degree of symmetry of FRW in order to 
write the curvature perturbations in such a way that they are still covariant tensors. We thus 
defin^ Rpupa as 



R (0) = 2 

A L pupa " 



H 2 + - 2 

a z 



(H 2 + H)hp a n y n p + perm. 



;ii8) 



Analogously, we define K^D = Hh^, which is relevant for what follows. Note that we have 
defined a background quantity by using the "full" tensor g pv and vector n p . From flllSp . 



4° ) = 2 (# 2 + ^) v-(# 2 + ^)(3 



R<® = Q{2H 2 + H+^- 

a 2 



(119) 
(120) 



11 Note that this is not exactly eq. (77) of Ref. (THJ: we have corrected few typos and chosen to use 



instead of <5° in the definition. 
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follow. Of the three dimension-four operators R 2 , R^R^ and R llU p cr R' J ' upcr , the latter can be 
substituted in favor of the other two and of the Gauss-Bonnet topologically invariant. 

G = R 2 - AR^R"" + R» vp „R? vp ° . (121) 

Another choice is to use quadratic combinations of the Weyl tensor C^ upa [3S] defined in 
eq. (1551) . which in our formalism is particularly handy because its background part vanishes 
on a FRW background. 

As R 2 can be expanded straightforwardly, let's consider R^ v R py in some detail. Follow- 
ing [18] we write 

R^BT = SR^SRr + 2Rf u R» v - R$R {0) . (122) 

The last term gives 

R (o) R (o)^ = 12 L H 4 + jj2 + 3H 2k_ + 3R 2jj + \ ^ (123) 

v a (X (X I 

and contributes to A(£). For the second term on the right-hand side of f)122p we find 

rWrT = 2(A-h) n^RT + (H 2 + H)R . (124) 



In turn (see e.g. 

n^BT = K 2 — K^K^ - V M (n"V F /) + V>"V/) , (125) 



The last two terms are dealt with in App. B of [15] : the first of the two gives something 
proportional to K\J— g 00 , which can be expanded in powers of 5g 00 . The last vanishes in 
unitary gauge. In summary, we find that an operator of the form A((j))R lJLV R tlv — > A{t)R lul R ,lv 
contributes to the terms in fl3J as follows (to simplify the formulas we set k = in the 
following) : 



(126) 



fit) = 2A(t) (3H 2 + 

A(t) = AA(t) (9H 4 + 5H 2 + 3H 2 H^j + 2H(AH) ' + 2{AH) 

c(t) = -UH(AH)- - 8AH 2 - 2{AH)" , 
M 2 4 = UH(AH)- + 8 AH 2 + 2{AH)" , 
fh% = -(AH)', M 2 = -Ml = 2{AH) , \ 2 = 2A. 

B Newtonian gauge 

After re-introducing n with the Stueckeberg trick, we are free again to choose the most ap- 
propriate gauge. The Newtonian gauge is particularly suited for late-time cosmology. In this 
gauge the metric is in the form (we assume hereafter a spatially flat FRW background) 

ds 2 = -e 2 *dt 2 + a 2 e- 2 ^5 ij dx i dx j . (127) 

It is not too long to calculate exactly the Christoffel symbols, 

r° = $, r^ = ^$, I% = a\e-™- 2 *{H-V), (128) 

r* = a~ 2 e 2 * +2 * d& , = SijiH - , T) k = 5 jk d& - 5^ - 5 ik d^ . (129) 
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With the aid of the above expressions we can calculate the extrinsic curvature defined 
in ([S]) as K^y = h°V c n u . By (JSJ), we have that in Newtonian gauge n p = — ^e*. Therefore, 



V p n u = -8 M e*d p ® + r°„„e* . 



(130) 



On the other hand, the projector h pu gives, for the space-space components hy = 5ija 2 e~ 2 ^ 
and zero otherwise. By using the expressions of the Christoffel symbols, eqs. fll28l) - fll29p . we 
finally get that the only non zero components of K pv are 



We also get 



K i:j = a 2 8ije- 2 *-*{H - V) ~ a 2 'S^H -V-H®- . 
5K = -3(tf + HQ) . 



(131) 
(132) 



B.l Gravitational action 

The Ricci scalar can also be written exactly in a rather compact form: 

R = 6 e~ 2 *(2H 2 + H -H$- AH^> + ^ + 2^ 2 - * ) 

- 2(T 2 e 2 * [V 2 $ - 2V 2 * + (V$) 2 - V$W + (V*) 2 ] 



(133) 



where V is a 3-d gradient and V 2 the Laplacian. 

It is useful at this point to write the Einstein-Hilbert term at quadratic order. After few 
integrations by parts we obtain 



-gR 



quadratic 



2(W) 2 -4V$W 



6a^ 



3,T r 2 



3a 3 



AHm + H 2 {<$> + 3*) 2 - 6* 2 (3# 2 + H) 



Finally we can multiply the Einstein-Hilbert action for the f(t) coefficient as in 
Stueckelberg trick f(t) produces a ir field, which therefore mixes with gravity, 



-gf(t)R -> 



/(*) + t(t> + 2^ 



R 



quadratic 



af 2(V*) 2 - 4V$V^ + 2afn (2V 2 ^ - V 2 $) - 6a 3 (/* 2 + fn^f) 



- 3a 3 / ^2$^ [2H + f/fj + ($ + 3^) 2 (# 2 + #/// 
•M' 2 ((i// 2 • HI ■ •'//./• / + /// - ,/'"". / 2 ) 

- 6a 3 /vr \H^> + 4#* + (2H 2 + + 3*)1 + a 3 fn 2 (QH 2 + 3if) . 



(134) 

By the 

(135) 

(136) 

(137) 
(138) 
(139) 



The first line, f)136p . is the most important because it contains the kinetic terms. This is the 
piece of the action that survives in the Minkowski limit^E The latter is obtained by posing 
H — H — f — 0. The reason why the /-terms in line f 11 3 6 j) survive in the Minkowski limit 
is that kinetic term for ir is multiplied by the coefficient c (see e.g. eq (]68p ). which means 



12 In Einstein frame this would be called the "decoupling limit" . In the Jordan frame however, as apparent, 
7r and gravity are kinetically mixed and thus not decoupled by taking the high-energy limit. 
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that the canonically normalized field is n c ~ irc 1 ^ 2 . In taking the Minkowski limit, the ratio 
/ '/c 1 / 2 should be kept constant and therefore also the combination fir. More explicitly, the 
Minkowski limit reads 

J V^gf(t)R ki = ic J 2/(W) 2 -4/V$W + 2/7r(2V 2 ^- V 2 $) -6(f^ 2 + fn^). (140) 
B.2 Linearized Einstein equations 

We want to derive the linearized Einstein equations in Newtonian gauge f)127p for action (|66|) . 
In this gauge the linearized Einstein tensor components read 

Goo = 3# 2 - QH^ + 2cT 2 V 2 * , 

= - a 2 [3H 2 + 2H]5ij (141) 
+ [2H(Q + 3*) + 2^ + 2(Q + ^)(3H 2 + 2H)]5 ij + V 2 ($ - ^)c% + didj^f - Q) , 
G 0i = 2%$ + HQ). 

Of course, since we want to study perturbations in Newtonian gauge, we have to reintroduce 
the scalar degree of freedom n via the Stueckelberg trick. It is handy to do so at the level of 
the action f lBB]) . The contribution to the Einstein equations of the first three terms is given 
in eq. (I14p . On the other hand, the operator M| affects only the (00) component yielding 
4M 2 4 ($-7r). 

In order to find the contribution from the operator fn\ we first reintroduce the 7r in the 
action, 

Sg 00 SK ->■ (g 00 + 1 + 2g^d ll ix){K - 3H + 3Hn - a" 2 V 2 vr) , (142) 
and we write the variation of K as 

5K = 5(V a n a ) = V a (n p 8g^) + I-^V^n^ 00 ) - \g a ^V a {5g^) . (143) 

Then, upon integration by parts, we obtain 



1, 



-y a (ml(g 00 + 1 - 2tt)) (n a n^n u + 2nJ™ - n a g^) , (144) 

where, at first order in metric perturbations, K — 3H = — 3( 1 f + HQ), g 00 + 1 = 2$ and the 
n a 's appearing explicitly only need to be computed at zero order in the perturbations. 

Taking into account these last two terms that we have discussed and assuming m\ constant, 
the linearized Einstein equations read 

(00) component: 

2/(a" 2 V 2 * - 3HV) + 3f(-Q + H 2 n + Hn - a~ 2 V 2 7r/3) + 3fHn 



M 2 



2c7r-(c + A)7r-2A$ + 4M 2 4 ($-7r) (145) 



3(tf + HQ) - 3H-K + a- 2 V 2 vr + 3H(Q - tt 



8T 00 . 
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By using the equations of motions fl24l) . f[25]) . (1211) and 5T 00 = Sp m + 2p m Q we get 



Mi 



2/(a~ 2 V 2 ^ - 3H(4> + H$)) - f(3(V + #$) + 3#7r + a _2 V 2 7T + 2#($ - tt)) 



/($ - tt) + (p D + p D )(Q - tV + 3#7r) + 4M 2 4 ($ - tt) 



3(tf + if$) - 3Hn + a- 2 V 2 vr + 3#($ - tV 



(146) 



(ij)-trace components: 



Ml 2f 



y + H<5> + 3H^ + (3F 2 + 2H)($ + + V 2 ($ - ^)/(3a 2 

+/ 6 + 2* + 4#($ + - (3# 2 + 2#> + 2V 2 7r/(3a 2 )j - 2#(/tt)' - (/tt 
+ 2c($ - tt) - 2^(A - c) + (A - c)tt - m?[$ - tt + 3#($ - tt)] = oT fe fc /(3a 2 

Again, by use of the equations of motions, 



Ml 2f 



ty + H<l> + 3H^ + (3# 2 + 2H)$ + V 2 ($ - ) /(3a 2 
+/ 2(# + HQ) + ($-#) + 3#($ - tt) - 3# 2 tt + 2V 2 7r/(3a 2 
-Pdtt + (pd +Pd)($ - tt) -m?[$ - tt + 3#($ - tt)] = 5p m . 

(ij)-traceless components: 



2/($ + vD)} 
(147) 



(148) 



Ml(d i d j --5 ij V 2 ) [/(*-$)-/ 



TT 



(149) 



(Oi) components: 

2/(* + + (/ - + /($ - 7r)l -2(p iJ +p JJ )9 i 7r-2m?9 i ($-7r) = 5T l0 . (150) 



M 2 «9, 



C Field redefinitions of the metric 

C.l Perturbative field redefinition in the Einstein frame 

This subsection of App. Oand Sec. 15. II are the only places in the paper where Einstein-frame 
quantities are not denoted with a hat nor with the subscript E and Jordan-frame quantities 
explicitly carry a subscript J. Here we show that the action 



S = J d A xy^g 
can be rewritten as 
S = d 4 Xy/^g 



M?R 



c(t)g 00 - A(t)+F^(Sg m ,6K. 



\wy • • • 1 



t) 



+ S m [gj V ,A 



(151) 



M 2 R 



~c(t)g 00 -A(t) + F^(5g 00 ,5K, u ,...;t) 



+ S m \f{t)g^M ■ (152) 
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Above, gf is defined as in eq. ( 11001) . i.e. 



9j 



/(*)<r+<r 



Sg°^g 0u [io(t) + n(t)6g 00 + ...} 



+ ±m^{t)g^5g™K + Knf{t)8K^8K; + ... 



(153) 



Consider now the field-redefinition of the metric 

_). = + f^SF^. (154) 

In terms of the new field g^, the matter action has already the form that we want, 

S m [gTM = S m [f{t)r v M ■ (155) 

Let us now concentrate on the DE part of the action and in particular on the first three terms 
inside the square brackets in eq. (11511) . To see the effects of a perturbative field redefinition 
consider the linear term in the functional Taylor expansion, 



Slg^-f-^SF^} ~ S[g^] 



5S 



-r\t)T(t)h(t)g 



~oo 



(156) 



where the ellipses denote second-order terms or higher. For the second equality we have used 
the background Einstein's equation, eq. (I103p . and for 5F^ V we have replace the second term 
of the Jordan frame metric, eq. (I153p . Thus, this Lagrangian transforms in itself plus a linear 
term in g 00 proportional to the parameter /?i(t), which redefines the function c(t), 



c(t) c(t) = c(t) + \f{t)j^ 



(157) 



plus second-order or higher terms which do not involve the matter energy-momentum tensor. 
The higher-order terms F^ in eq. (11511) are clearly affected by the field redefinition. However, 
their structure remains unaltered, as they can always be rewritten as polynomials of the 
geometrical invariant objects 5g 00 jSK^, etc., following [17]. In conclusion, the action after 
the field redefinition is eq. (I152p . 



C.2 From Jordan to Einstein 

In this appendix Einstein frame quantities are denoted with an index E or with a hat. In order 
to study the stability of our DE theory we consider action (j4j) keeping terms up to second 
order in the perturbations, 



5 



DE 



d A Xy/—g 



M 2 

-ff(t)R-c(t)g™-k(t) 



ML, 

+^r($g 



00\2 



2 2 2 



(158) 
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In order to switch to Einstein frame it will be useful to review some standard "vocabulary" of 
conformal transformations, 



9V 
^9 



f 1 9 l iu 

2 r 



f 



9 , 



R = f ( R E + 3D In / - -g^ In fd v In / 



(159) 
(160) 

(161) 



from which we get 



» 4 v /Z^oo / 



fl/UxT> 3 5g™ If 



(162) 



where h is the metric projected on t = const, surfaces. The problem with g 00 so defined is that 
its background value is — f~ l instead of —1. This is related to the fact that by conformally 
transforming the metric components (eq. (I159|) ) without changing coordinates one looses the 
gauge choice e/oo = g 00 — — 1. While moving to the Einstein frame it is thus useful to redefine 
the time coordinate (He — f 1 ^ 2 dt in order to have the metric in the usual form. This was done 
systematically in Sec. 15.21 After such a time redefinition we find 



5 



DE 



d 4 XE 



M 2 

—^-Ke — ce9 — A E 



A/f 4 777 3 A m 2 -s M 2 

iW 2 ^~00\2 m U;00tr> M 2 fr >2 IV1 3 



(163) 



with 



/ 2 M 2 4 



£ 2 

(•3/2^3 



f'c E = c+-Mi 



f 



f 2 A E = A 



Mi 



l^(f/f)-^M 2 (f/f) 2 -^Mi(f/ff 



f 2 nf\ = m\ + \m\ (f/f) + Hi 2 , (f/f) , fM, 



mi, 



fM. 



Mi, 



(164) 
(165) 
(166) 



where, on the right-hand side, a dot denotes the derivative with respect to J-frame time. 
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